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Phase diaeram of 5 = I XX Z chain with NNN interaction 
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Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 
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We study the ground state properties of one-dimensional XXZ model with next-nearest neighbor 
coupling a and anisotropy A. We find the direct transition between the ferromagnetic phase and 
the spontaneously dimerized phase. This is surprising, because the ferromagnetic phase is classical, 
whereas the dimer phase is a purely quantum and nonmagnetic phase. We also discuss the effect of 
bond alternation which arises in realistic systems due to lattice distortion. Our results mean that 
the direct transition between the ferromagnetic and spin-Peierls phase occur. 
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Quantum spin chains have attracted a considerable 
amount of attention over the past decades. The main 
reason for it is the dominant role played by quantum 
fluctuation and frustration in these systems One of 
the simplest model of them is the spin-i XXZ chain 
with next-nearest-neighbor interaction, described by the 
Hamiltonian 



H = J '^{hi^i+i + ah 



— Oi Jj 



'i,i+2) 



SfS^ + ASfS^ 



(1) 



In the following, we set NN coupling J > 0. This system 
can also be considered as a coupled chain with zigzag in- 
terchain interaction. Recently, the physical interest for 
this model has been renewed, with relation to experi- 
ments; it is discussed that this system describes a spin- 
Peierls compound CuGeOa [||. 

The ground state properties are known for several 
cases. When a = 0, it is integrable by the Bethe ansatz 
For A < — 1, the system becomes the ferromagnetic 
phase, while for A > 1, it becomes the Neel phase. In 
the case of XY-type anisotropy (—1 < A < 1), quan- 
tum fluctuation destroys the long range order even at 
the zero temperature. The ground state is characterized 
by a gapless excitation and algebraic decay of spin corre- 
lations (spin fluid phase). On the a = ^ line, the other 
exact results are obtained. At isotropic point (A = 1), 
ground state is purely dimerized |^] and there exists a 
finite energy gap above the ground state Q. Recently, 
it is proven that the ground state is still dimerized for 
A > — i, while for A < — i a ferromagnetic phase ap- 
pears g. 

In this letter, we determine the phase diagram of this 
system in the region < a < ^, A<0 (see FIG. |l|) by 
use of diagonalization of finite chain Hamiltonian. This 
model has been investigated by many authors ||^,^ . The 
phase diagram of the other region A>0,0<a<i 
has been obtained by Nomura and Okamoto From 
FIG. |l|, we see that the transition between the ferromag- 
netic phase and the dimer phase occurs in rather broad 



region in the phase diagram. This result is surprising, 
because the ferromagnetic state is classical, whereas the 
dimer state has a purely quantum nature. The magnetic 
properties may drastically change at the ferro-dimer tran- 
sition line. 

First, we explain how to determine the boundary of fer- 
romagnetic phase. The transition from the = ±L/2 
(Ferro) state to the — (dimer or spin fluid) state 
occurs via a level crossing. Thus, we obtain the bound- 
ary of ferromagnetic phase by the level crossing between 
the fully magnetized state {S^ = L/2), which have the 
energy eigenvalue E — -^(1 + a)L, and the singlet state 
(5^ = 0) Then, we compare the numerical result 
with the variational one. At a = ^,A > — i, ground 
state is purely dimerized state 

|$i) = [l,2]...[L-l,L], |$2) = [2,3]--.[L,1], (2) 

where is a singlet pair of spin i and j. Therefore, 
near a = i, we take a trial wave function as |$«Qr) = 
|$i) -f /3|$2), and estimate the ground state energy Eyar 
variationally. We find 
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FIG. 1. Phase diagram of the S — ^ 
with NNN interaction. Hatched area denotes gapless spin 
fiuid phase. The dimer-spin fiuid boundary is denoted as 
K = 1. Open symbols are related to the instability of 
spin fluid phase; K = 2 (o) and K — 4 {o). The lines 
K — 1, 2, 4 terminate at the boundary of ferromagnetic phase, 
(a. A) = (0.393, -0.545), (0.296, -0.608), (0.214, -0.683), re- 
spectively. 
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^(A + 2)(i-.), (3) 



and space inversion is expressed in a similar way, 



where A and B are positive coefficients of order 2^^/^. 
Tire ferro-dimer phase boundary is estimated from 

Eferro = Eyar, which yiclds, for L —^ oo, 



2ac + 3 



It gives, for — > ^, 



1 1 

-2 + 4'"^^ 




(4) 



(5) 



It coincides with our data up to the first order of ac — ^ • 
The dimer-spin fiuid transition is known to be of 
Berezinskii-Kosterlitz-Thouless (BKT) type It has 
been difficult to determine the BKT point and critical 
indices numerically, since transition occurs with essen- 
tial singularity |pl]| . In addition, there appears logarith- 
mic correction as a finite size effect which converges very 
slowly. These features make it difficult to determine the 
transition point accurately. In fact, a naive application 
of finite size scaling methods often lead to a misleading 
conclusion jl|] . Recently, an efficient method, level spec- 
troscopy, has been proposed to resolve these difficulties, 
based on the renormalization group argument and sym- 
metry consideration ^ . 

Following the standard steps of Jordan- Wigner trans- 
formation and bosonization, one obtains the quantum 
sine-Gordon Hamiltonian, which describes the large dis- 
tance behavior of the Hamiltonian (|l|) [|l^,^ , 



H=^l dx [vKiTTUf + ^{d<l>r} 

dx cos VSip, 



'27ra2 



(6) 



where a is a short-distance cut off or a lattice spacing, v 
is spin wave velocity. 11 is momentum density field con- 
jugate to (p, [(j){x),Il{x')] — iS{x — x'). We introduce the 
field 9 dual to cj), defined as dxO — ttII. We compactified 
those fields a.s cj) = (l)+-^,0 = 6+-^. The spin operators 
are represented by 4>,9, 



S: 



X I a 



dx 



-6"="/'^ cos\/20(a 



(7) 



and is expressed as g^'^^. The symmetry operations 
are interpreted as follow: translation by one site 



Tr 

spin reversal 
T:< 



TT 

71' 



TT 



TT 

71' 



TT 

71' 



(8) 



(9) 



V2' 



(10) 



When Hff, — 0, i.e., the Gaussian, all correlation func- 
tions can be calculated. The excitation spectrum is gap- 
less (massless) and characterized by critical behavior. 
When ^ 0, we follow the renormalization group ar- 
gument. For K ~ 1 and ~ 0, we obtain following RG 
equations |14|, 



de 
dym 

d£ 



(11) 



where yo{i) and y^{i) are related to bare couplings as 
K = 1 + ^yo{0), y^{0) = y^, and the scaling factor £ is 
related to the system size L as £ = InL. For j/q > \y<p\i 
y^ term is irrelevant and it is renormalized to zero for 
L — > oo. Then, ground state properties are the same 
as the case of = except K is renormalized, which 
correspond to spin fluid phase in the spin system. While 
for yo < ly^li ?/0 term becomes relevant and it grows large 
at large distance limit. Then, quantum fluctuation is 
suppressed and spontaneous symmetry breaking occurs. 
The system becomes a Neel or dimer phase depending 
on the sign of the bare coupling constant y^ The 
marginal case, yg = \y^\, is just critical line and both 
couplings are renormalized to zero. Thus, the critical 
point is identified as a point at which the coupling K is 
renormalized to one. 

In the spin fluid region, large distance physics is de- 
scribed by the Gaussian model with (renormalized) cou- 
pling K. Then, the system possesses conformal symme- 
try [|l5|. As a consequence, there exists a state, in a 
periodic chain of length L, associated with each primary 
operator of scaling dimension x and spin s, whose exci- 
tation energy and wave number is given by p6[ , 



. „ 27rw , 27r 

A£; = X, k = — s. 

L L 



(12) 



In our case, the vertex operators [ pT[ are most significant, 

(13) 



n.ra = 6 C , 



which have, 



1 / 2 m2\ 

Xn,ni = -\Kn + I : Sn,™ = "^'^^ (14) 

where n, m are integers. 

From ^ and (|l|), we see that, at critical point 
K — 1^ some excitations become degenerate, like On eel = 

Olfl + O-1,0j Odimer = Oifi — 0-1,0, ^nd Odoublet = 

Oo,±i; these excitations have scaling dimension x — ^. 
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In the spin system, corresponding excitations are classi- 
fied as, = 0, k = TT, P = -1 (Neel), S'f, = 0, k = 
TT, P = 1 (dimer), 5f = ±1, fc = tt, P ^ -1 (doublet) 
[cf. (|^)-([lO|)]. Thus we can determine the transition 
point from the crossing of those excitations. 

In the sine-Gordon model at criticality, there appears, 
as a finite size effect, logarithmic correction of 0(1/ In L) 
from the marginally irrelevant field cos \/S(t' i thus 

the degeneracy of those fields splits. However, there re- 
mains some degeneracy. The splitting of four x = ^ fields 
forms triplet and singlet reflecting implicit SU{2) sym- 
metry on the BKT transition line |Q. Therefore, we can 
eliminate the logarithmic correction by choosing member 
of triplet, as O dimer and O doublet, and we can determine 
the BKT point accurately. Besides the logarithmic cor- 
rection, there remains a correction of 0(L^^), originated 
from the irrelevant fields (x = 4) which is not included 
in Hamiltonian (^ [p^ . 

As for the critical indices, logarithmic correction can 
be eliminated in the averaged scaling dimension such as 
{xNeei + Xdimer + '^x doublet) / sincc the ratio of lead- 
ing logarithmic correction is known as C^eei '■ Cdimer ■ 
Cdoubiet = 3 : — 1 : —1 reflecting the implicit SU{2) sym- 
metry of the BKT transition. 

Figure |^(a) shows a size dependence of crossing points 
ac{L). We see that the leading correction in adL) is 
0{L~^) as expected. We check the averaged scaling di- 
mension mentioned above, shown in FIG. ^(b), which 
converges to the value ^. We also estimate the confor- 
mal anomaly number c along the BKT line, from the 
finite size correction to the ground state energy pT] , 



Eo{L) 



TTV 



(15) 



Our estimate is, c = 0.9956, 0.9958, 0.9960 for A = 
0.0, —0.2, —0.4, respectively. We conclude that, together 
with the averaged scaling dimension, the system belongs 
to the same universality class as sine-Gordon model. 

In realistic systems, we should consider the lattice dis- 
tortion or the staggered field. The spin fluid phase be- 
comes unstable against such perturbations, and the sys- 
tem is no longer gapless for i^T < 4. The line K = 2 is 
related with the spin-Peierls instability. 

For instance, we consider the bond alternation 



Halt.=Sj2i-^ySr Sj + U 



which has a following bosonized form 



25 



(27ra)2 



dx sin 



(16) 



(17) 



where g ^ 6 + 0{S^), because it must be invariant under 
the change 6^—5 followed by one site translation (f> 
(/) -|- It explicitly breaks the symmetry of one site 
translation, though the system remains invariant under 



the two site translation. One can see that this form is 
also invariant under the spin reversal T. The operator 
sin\/2(/) has scaling dimension x = K/2, thus it becomes 
relevant and causes a gap ior K < 4. 

Let us explain how the K = 2 line is related to the 
spin-Peierls instability. When K < 4, S = is just a 
critical point, and correlation length ^ diverges as 



which is followed by the RG equation for S 
dS{i) 



d£ 



{2~x)Sii)+OiS^). 



(18) 



(19) 



To see how the ground state energy changes with the 
strength of bond alternation, we consider the following 
susceptibility 



X 



(20) 



where Eq is the ground state energy. For (5^1, we flnd 

x/L-e-'^'-S^-'. (21) 

Therefore, the ground state energy per site eg behaves as 

eo(<5)«eo(0)-A.<54/(4-^), (22) 

where A is some positive constant . The ground state 
energy decreases as 5", a = 4/(4 — K) due to the for- 
mation of singlet pairs. When i^T < 2, it decreases the 
ground state energy by a sufficient amount to compensate 
the energy cost of lattice distortion 5^, therefore spin- 
Peierls instability occurs. 
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FIG. 2. (a) Size dependence of crossing point ac(A) defined 
as A_Edou6!et(Qc, A) = AEdimerictc, A) and (b) the averaged 
scaling dimension (xjveei + ^ dimer + 'iXdouUet)/4 obtained at 
BKT point [K = 1), for A = 0. 
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To obtain the K = 2 and K = A line, we impose the 
twisted boundary condition (TBC) with boundary angle 
TT which effectively shifts the value n by an amount ^, i.e., 
n changes to n+ i under twisted boundary condition [ p3[ . 
The excitation AE'^q'^ crosses with AiJo.i or Ai?o,2 at 
K — 2 or K ^ 4, respectively. We also calculate the ratio 



Ai?i,o/A£;o,i = a::i,o/a;o,i 



(23) 



under periodic boundary condition. According to (14), 
it should be equal to K^. It can be used to check the 
result obtained by the method mentioned above, and we 
confirmed within 1% accuracy (see TABLE. |). 

The K — 1, 2, 4 lines terminate at ferro-spin fluid 
transition line and they do not merge at the same point. 
It indicates that the value K is finite at ferro-spin fluid 
transition line except at a = 0. This result is somewhat 
unexpected. In the spin fluid region, spin correlation 
behaves as 



K 



+ const. X (-1)^/'^ 



~K 



(24) 



At first sight, appearance of ferromagnetic order may be 
interpreted as K oo where the uniform part of spin 
correlation diverges and the staggered one vanishes [Q. 
But our result suggests that this picture is incorrect for 
a 7^ 0, i.e., K remains finite at ferro-spin fluid phase 
boundary except ata = 0, A = — 1. 

In summary, we determined the phase diagram of 
one-dimensional XXZ model with next-nearest neigh- 
bor coupling a and anisotropy A in the region Q < a < 
i, A < (FIG. d). A striking feature of the phase dia- 
gram is that the ferromagnetic phase and the dimer phase 
share a boundary of finite length in the phase diagram. 
In a realistic situation, where we should include the lat- 
tice distortion in addition to the spin-spin interaction, 
we will observe the transition between the ferromagnetic 
order and the spin-Peierls order, which can be observed 
in a broader region (K < 2 in FIG. I). 

The numerical calculation in this work was based on 
the computer code TITPACK ver. 2 developed by Pro- 
fessor H. Nishimori. 



TABLE I. Spin wave velocity v, the conformal anomaly 
number c, and the ratio 3;i,o/a;o.i on K = 2, 4 lines. 



Extrapolated value 



(A, a) 


V (Exact'') 


c 




[20 


( 0.00,0.00 ) 


1.000 (1) 


1.000 


4.00 




(-0.20,0.097) 


0.7463 


1.000 


4.00 




(-0.40,0.193) 


0.4918 


1.000 


4.00 


[21 


(-0.60,0.292) 


0.2402 


0.996 


4.01 


[22 


(-0.707,0.00) 


0.4731 (0.4714) 


1.000 


16.0 


(-0.685,0.10) 


0.3272 


0.997 


16.0 




(-0.682,0.20) 


0.1990 


0.984 


15.9 


[23 


^ Exact value of v 


is known for a = 0. 
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